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ABSTRACT 

Motivated by the observed correlation between black hole masses Mbh and the veloc- 
ity dispersion a of host galaxies, we develop a theoretical model of black hole formation 
in galactic bulges (this paper generalizes an earlier ApJ Letter). The model assumes an 
initial state specified by a a uniform rotation rate Vt and a density distribution of the 
form p = alg/27rGr'^ (so that Oofr is an effective transport speed). The black hole mass 
is determined when the centrifugal radius of the collapse flow exceeds the capture radius 
of the central black hole (for Schwarzschild geometry). This model reproduces the ob- 
served correlation between the estimated black hole masses and the velocity dispersions 
of galactic bulges, i.e., Mbh ~ lO®M0((T/2OOkms-i)^, where a = V^a^g. To obtain 
this normalization, the rotation rate 17 ~ 2 x 10^^ rad/s. The model also defines a bulge 
mass scale Mb- If we identify the scale Mb with the bulge mass, the model determines 
the ratio fiB of black hole mass to the host mass: fiB ~ 0.0024 ((T/200kms~^), again in 
reasonable agreement with observed values. In this scenario, supermassive black holes 
form quickly (in ~ 10^ yr) and are born rapidly rotating (with a/M ~ 0.9). This paper 
also shows how these results depend on the assumed initial conditions; the most impor- 
tant quantity is the initial distribution of specific angular momentum in the pre-collapse 
state. 

Subject headings: black hole physics - galaxies: nuclei - galaxies: kinematics and dy- 
namics 



1. INTRODUCTION 



During the past decade, the observational evidence for massive black holes has crossed a 
threshold of reliability and black holes are now considered to be discovered. Almost every galaxy 
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is thought to harbor a supermassive black hole anchoring its center (e.g., Magorrian et al. 1998; 
Kormendy & Richstone 1995). Our own Milky Way galaxy contains a modest central black hole, 
with its mass estimated at M^w ~ 3 x 10'^ Mq (e.g., Genzel et al. 1996, Ghcz et al. 1998). The 
properties of these black holes and their connections to their galactic hosts are currently the subject 
of intensive investigation. 

A striking aspect of the black hole-galaxy connection has been recently reported. Two com- 
peting groups have observed a relationship between the velocity dispersion a of the host galaxy and 
the mass Mbh of its central (supermassive) black hole (Gebhardt et al. 2000; Ferrarese Sz Merritt 
2000). This correlation can be written in the form 

Mbh = Mo(a/200 kms-^)T, (1) 

where the two observational teams find slightly different values for the constants in this scaling 
relation. The exact values derived from the data depend on the fitting procedure and are sensitive 
to the exclusion of outlying points (see Merritt &: Ferrarese 2000, 2001). A recent in-depth analysis 
(Tremaine et al. 2002) finds that 7 = 4.02 ± 0.32 with a mass scale Mq = 1.3 x 10^ M©. In any case, 
the observed correlation is remarkably tight: The observed scatter in black hole mass M^h at fixed 
dispersion a is less than 0.30 dex (about a factor of 2). Furthermore, the observed relation appears 
to be independent of the Hubble type, profile type, or galactic environment. Previous observational 
surveys have found correlations between the black hole mass and bulge luminosity (Richstone et al. 
1998; Magorrian et al. 1998; van der Marel 2000; Kormendy & Richstone 1995; see also CaroUo, 
Stiavelli, &; Mark 1998), but the relation [1] appears to be far more robust. This observed scaling 
relationship provides an important constraint on theories of galaxy formation and bulge formation. 
Such theories must ultimately account for the production of supermassive black holes at galactic 
centers and the observed relationship between black hole mass and galactic velocity dispersion. 

In a previous paper (Adams, Graff, &; Richstone 2001; hereafter Paper I), we presented a 

theoretical model for black hole formation during the collapse and formation of galactic bulges. 
This model uses an idealized treatment to describe the collapse of the inner part of protogalaxies. 
The initial state is assumed to have a density distribution of the form p oc (like that of a 
singular isothermal sphere) and a uniform rotation rate; the initial conditions are characterized by 
an effective transport speed Oeff and the rotation rate J7. As the collapse develops, material falls 
inward from ever larger starting radii and carries larger amounts of specific angular momentum. 
The black hole mass is determined when the centrifugal radius of this collapse flow exceeds the 
capture radius of the black hole growing at the center. In spite of its idealized nature, this simple 
model correctly accounts for the observed Mbh — cr relation (equation [1]) with no free adjustable 
parameters and also predicts the observed ratios of black hole mass to the host (bulge) mass scale. 
Because of this preliminary success, the model deserves further exploration, which is the subject of 
this present work. 

We note that many other theoretical models have been developed to explain the observed rela- 
tionship between hole mass and galactic velocity dispersion (e.g., see the recent review of Richstone 
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2002). A semi-analytic model of merger-driven starbursts with black hole accretion (Haehnelt & 
Kauffman 2000; Kauffman & Haehnelt 2000) provides a correlation of the observed form (with the 
proper choice of the free model parameters). Several models are based on the idea that black hole 
accretion can influence star formation and gas dynamics in the host galaxy; this feedback can occur 
through ionization, mechanical work, and heating (e.g., Ciotti k, Ostriker 1997, 2001; Blandford 
1999; Silk & Rees 1998). The model of Blandford (1999) predicts that Mbh < r]a^ whereas the 
model of Silk & Rees (1998) implies Mbh oc . The idea that the black hole mass is limited by disk 
accretion has been explored by Burkert & Silk (2001). Before the observational correlation was 
discovered, Daniel & Loeb (1995) argued that the seeds for quasar black holes could originate from 
the collapse of low angular momentum regions. Finally, the accretion of collisional dark matter 
indicates a scaling relation of the form oc a^~^-^ (Ostriker 2000). 

This paper is organized as follows. In §2, we review and extend the model presented in Paper 
I, and describe the orbital infall solutions in greater detail; we also generalize the model to include 
continued infall onto the black hole at late times. In §3, we consider mass accumulation onto 
the black hole through disk accretion, more general initial conditions, the effects of mergers, and 
nonzero quadrapole moments in the initial conditions. We conclude, in §4, with a summary and 
discussion of our results. 



2. THE ROTATING COLLAPSE MODEL 

In this section, we review and expand upon the basic model of black hole formation during 
the collapse of a forming galaxy (see Paper I). In this context, we examine the collapse of the inner 
part of a region that will ultimately form the bulge of a galaxy. 



2.1. Initial Conditions 

The calculation starts at the time of maximum expansion for the main body of the bulge. The 
main characteristics of the model can be summarized as follows: 

[1] All of the matter participating in the collapse - including baryons, dark matter, and stars 
- are unsegregated. In particular, we assume that all of the collapsing matter has the same initial 
distribution of specific angular momentum, as this distribution is the key ingredient in producing 
black holes with the observed properties. However, this model does allow for the possibility that 
additional material, perhaps some of the dark matter, does not participate in the collapse (see 
below for further discussion). 

[2] The mass and density distributions in this region take the form of a singular isothermal 
sphere even though the system is not in virial equilibrium. More specifically, the initial density and 
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mass distributions are assumed to have the form 

^"^^^ = 2^ M(r) = -^r. (2) 

The effective transport speed that specifies the initial conditions is related to the isotropic 
velocity dispersion a according to cr = v^o-cfr, where a is the velocity dispersion of the final state. 
This relation results from converting half of the original potential energy into kinetic energy, with 
the overall radius of the structure shrinking by a factor of 2 (see also below). For gaseous material, 
the transport speed Oeg plays the role of the sound speed. For any dark matter participating in 
the collapse, the intrinsic velocity distribution of its initial state is highly ordered with width Sv <^ 
Oeff ~ cr (by assumption; see also below). 

[3] This region is slowly rotating like a solid body (e.g., due to tidal torques) at a well-defined 
initial angular frequency J7. In this starting configuration, both dark matter particles and parcels 
of baryons that are initially located at radius r^o have initial angular momentum j = ^{r^o sin^o)^; 
where is the (initial) polar angle in spherical coordinates. 

[4] The central region of the collapse flow successfully produces a "seed" black hole in the 
earliest phases of evolution. The mass of this starting black hole may be much smaller than the 
large (Mbh ~ W^Mq) black holes of the final states. Once a black hole has condensed out of the 
galactic center, it will grow according to the calculations of this model. The initial seed black hole 
could form by the collapse of the densest (central) part of the perturbation or could be primordial. 

The initial state is characterized by two physical variables: Oeg and Q. We stress that these 
quantities are not free adjustable parameters, but rather can be specified - or at least constrained 
- by observations. First, we note that the initial transport speed OeS is directly related to (but not 
equal to) the final velocity dispersion of the final system. In collapsing regions with no dissipation, 
the virial theorem implies that the scale length of the mass distribution drops by a factor of 2 from 
the point of maximum expansion (see Binney & Tremainc 1987). Observational considerations 
(e.g., the "flat rotation curve conspiracy") suggest that dissipation does not greatly alter the final 
value of the dispersion a. This argument implies that the observed velocity dispersion a is related 
to the initial transport speed UeS of the protogalactic material through the relation = 2agg. 

Throughout this paper, we adopt a fiducial value for the starting rotation rate J7 = 6 x 
10~^Myr~^ = 2 X 10~^^ rad s~^, which ultimately provides the observed normalization for the 
-^bh — c relation. The scatter about this fiducial value will produce a corresponding scatter in 
the theoretical Mbh ~ ^ correlation. Although the initial rotation rates of the inner portions of 
galaxies in their pre-coUapse states are no longer observable, this adopted value is in reasonable 
agreement with expectations. For a ballpark estimate, we can consider the fundamental plane 
(Binney Sz Merrifield 1998), which provides a relationship between the half-light radii of galactic 
bulges and the corresponding velocity dispersions. For a typical value of the velocity dispersion 
a = 200 km s~^, the effective radius Re of a galaxy on the fundamental plane is about 3.5 kpc. 
In order of magnitude, we expect that ^ a /Re ~ 2 x 10~^^ rad/s. Rotation rates of this 
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order of magnitude are also consistent with those expected from observed rotational velocities in 
galactic bulges (sec, e.g., Figure 4.6 of Binney &; Treamine 1987; Binney &; Merrifield 1998; Jarvis 
& Freeman 1985; Wyse k Gilmore 1992). 

It is useful to compare our assumptions to the typical value of the spin parameter A for 
protogalaxies predicted by numerical simulations (where the angular momentum originates from 
cosmological torques - see Peebles 1993). Following Bullock et al. (2001), we write the spin 
parameter in the (slightly non-standard) form 



V^RVcir 

where J/M is the specific angular momentum, R is the outer radius, and Vdr is the circular speed of 
the protobulge structure. For our assumed initial density distribution [2] with solid body rotation, 
the specific angular momentum j = J/M = (2/9)r2i?^ and the circular speed at R is given by Vdr 
= \/2aeff. Using these results in the definition [3], we find that A = 1/9 for our assumed initial 
condition. For comparison, the values of A that are predicted for dark matter halos by numerical 
studies of structure formation have a mean value A = 0.04 — 0.05 (Peebles 1993; Barnes & Efstathiou 
1987; Bullock ct al. 2001). Although these values are a factor of 2 smaller than used here, they are 
evaluated for halo mass scales that are thousands of times larger than the protobulge mass scales 
used here. Notice also that if the bulge forms via collapse with dissipation, then the bulge spin 
parameter can be larger than the halo spin parameter (e.g., see White 1996). 

We can think of the initial conditions (oes , fi) in two different ways: First, we can consider the 
effective transport speed Oefr and the rotation rate O as adjustable parameters that can be varied 
in order to explain four quantities: the velocity dispersion cr, the bulge size scale i?, the bulge 
mass scale Mb, and the central black hole mass Mbh- On the other hand, we can use a and R to 
specify the initial parameters a^s and Vl. In this latter case, we arc left with a theory containing 
no adjustable parameters, but the theory must still correctly account for the bulge mass scale Mb 
and the black hole mass Mt,h as a function of a. 

As we show below, all material with initial conditions given by equation [2] follows a ballistic 
trajectory as it falls toward the central black hole. This result holds for gas, stars, and dark 
matter. Gas naturally takes on a centrally concentrated distribution and approaches the form of 
equation [2] ; the gas density obtains this form for the limiting case of hydrostatic equilibrium with 
an isothermal equation of state. However, the dark matter is somewhat less likely to assume this 
same starting state. Unlike gas particles, individual dark matter particles can have high angular 
momentum even if they are part of a larger structure with low (or zero) angular momentum. In 
other words, it is possible for the dark matter to display the overall density distribution of equation 
[2] , and for the structure as a whole to rotate slowly, and still have the individual particles possess 
too much angular momentum to be captured. In order for dark matter to fall into the central black 
hole, the particles must be extremely cold (with internal velocity dispersion 5v <^ a^s ~ cr), and 
the large scale streaming velocities Vs in nonradial directions must also be small {Vs ^ aeS ^ a). 
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Observations have shown that the central black holes of galaxies are strongly correlated with 
the properties of galactic bulges, but are more weakly correlated with other galactic components 
such as the disk or the dark matter halo. For example, the galaxy M33 has no bulge (but has both 
a disk and a dark matter halo), but does not contain a black hole. Although galactic bulges may 
contain some dark matter, they are probably not dominated by dark matter (e.g., Gerhard et al. 
2001). The micro-lensing optical depth of our own Milky Way bulge is high, ~ 2 x 10~^, indicating 
that it contains many baryons in the form of stars, stellar remnants, and brown dwarfs (e.g., Evans 
&: Bclokurov 2002). In fact, the baryonic content of our bulge is so high that it is difficult to 
construct dynamical bulge models using all of the baryons; very little of the mass budget is left 
over for a significant dark matter component. Dark matter may not play a dominant role in our 
galactic bulge, and, by extension, may not dominate the determination of the observed M^h — cr 
relation. 



2.2. Classical Orbit Solutions 

As the collapse proceeds, particles in the initial distribution fall towards the galactic center. 
Because the dynamical time scales monotonically increase with radius, infalling shells of material 
do not cross. The mass contained inside a given spherical shell, which marks a particle's location, 
does not change as the particle falls inward and hence orbital energy is conserved. In the classical 
(nonrelativistic) regime, the orbital energy is given by 



^ In if GM 



In this problem, we consider orbits that begin their collapse trajectories at large radii and then fall 
a long way toward the center of the galaxy. As a result, we can idealize these trajectories as zero 
energy orbits. 

For a given gravitational potential, we find the orbital solutions for material falling towards 
the galactic center; the same orbits apply to stars, dark matter, and parcels of gas. In our initial 
calculation (Paper I), the inner solution is derived using the gravitational potential of a point 
source. This form is only used in the innermost regime of the collapse flow where the potential is 
dominated by the forming black hole. In other words, this orbital solution derived here is valid 
over the range of length scales 

i?5 < r < roo , (5) 
where Rs is the Schwarzschild radius of the black hole and is given by 

2GM 

i?5 = ^. (6) 

In general, the black hole produced through this collapse process will be rotating so that its event 
horizon is not completely specified by the Schwarzschild radius Rs- Notice that at late times, long 
after black hole formation is complete, dark matter and stars will miss the black hole and continue 
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to trace through their orbits back out toward large radii; this behavior leads to an extended mass 
distribution and the potential is no longer well described by a point potential. At these later times 
of evolution, our solution loses its validity. In the present context, however, we only use the solution 
during the early phases in which the gravitational potential of the central region is dominated by 
the black hole. In addition, we note that relativistic corrections become important as r — *■ Rs- 

Since this potential is spherically symmetric, angular momentum is conserved and the motion 
is confined to a plane described by the coordinates (r, (j)) ; the radius r is the same in both the 
plane and the original spherical coordinates. The angular coordinate in the plane is related to 
the angle in spherical coordinates by the relation cos^ = cos 9 / cos^Oj where is the angle of 
the asymptotically radial streamline (see below). For zero energy orbits, the equations of motion 
imply a cubic orbit solution, 

where joo is the specific angular momentum of particles currently arriving at the galactic center 
along the equatorial plane. The second equality serves to define the parameter C,. Here, the 
trajectory that is currently passing through the position (r, /x = cos 9) initially made the angle 
with respect to the rotation axis (where /xq = cos^o)- 

For a given angular momentum, we can use equation [4] to determine the pericenter p, which 
marks the distance of closest approach for a parabolic (zero energy) orbit. Our assumption of 
uniform initial rotation implies that j^o = ^i^'^ sin^ where is the starting radius of the 
material that is arriving at the center at a given time. The pericenter can be written in the form 

^ ^ (roosingo)^^^^ ^ {GMfn''sm^9o 

^ 2GM 2GM 25afjj ' ^ ^ 

In the final equality, we have used M = M{r) as a label for by inverting the mass distribution 
of the initial state (equation [2]) to find roo=GM/2a'^g. As in star formation theory (Shu et al. 
1987; Cassen and Moosman 1981), we define the centrifugal radius Rc of the flow according to 

n^G^M^ 

which represents the radius of a circular orbit with angular momentum j^o for incoming matter 
falling within the equatorial plane. For motion in the equatorial plane, this radius Rc is twice as 
large as the pericenter p for a parabolic orbit with the same angular momentum. 

Given the orbital solution (equation [7]), we can find the velocity fields for the collapse flow. 



Vr = - 



■^)">-C(l-.S)}'^ (10) 

>j 1/2 

-/^')C , (11) 
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\ 1/2 

GM^ 



^^=(^^J {i-^^i){i-i^r"'C''- (12) 

Since Mj ^iid /xq are related through the orbit equation [7], the velocity field is completely 
determined for any position (r, 9) . 

The density distribution of the infalling material can be obtained by applying conservation of 
mass along a streamtube (Terebey, Shu, & Cassen 1984; Chevalier 1983), i.e., 

p{r, 0) Vr sin 6 dO dip = sin dOg dipo, (13) 

47r 

where M is the total rate of mass flow (inward) through a spherical surface (e.g., Shu 1977). 
Combining the above equations, we can write the density profile of the incoming material in the 
form 

The properties of the collapsing structure determine the orbit equation [7], which in turn determines 
the form of diio/dji. With the radial velocity given by equation [10], the density field is thus 
completely specified in analytical form (implicitly). 

Notice that we have ignored gas pressure in the collapse solution. Dark matter and stars 
always exhibit pressure-free (ballistic) behavior and our approximations are automatically justified 
for these components. Even for gas, however, the collapse flow always approaches a pressure- 
free form in the inner region. This (somewhat remarkable) characteristic follows from considering 
the gaseous portion of the collapse flow to be a scaled-up version of the collapse flows that have 
been studied previously for star formation theories (Shu 1977; Terebey, Shu, &: Cassen 1984; see 
also Adams 2000). In this case, the collapse of the initial state (with density distribution [2]) 
proceeds from inside-out and gas parcels in the central portion of the flow always approach ballistic 
trajectories. 



2.3. The Mass Scale for Galactic Black Holes 

With the collapse solution in place, we can estimate the mass of forming black holes. The 
collapse flow defines a critical mass scale Mc, which will be roughly comparable to the final black 
hole mass. In the earliest stages of collapse, incoming material falls to small radii r < Rs, where 
Rs is the Schwarzschild radius of the forming black hole. The mass that determines Rs is the total 
mass M = Mt that has fallen thus far, i.e., we assume that the black hole mass M^^ = M in this 
early evolutionary stage. As the collapse proceeds, incoming material originates from ever larger 
radii and carries a commensurate increase in specific angular momentum. The centrifugal barrier 
of the collapse flow thus grows with time. 
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If the pericenter p is sufficiently small, ballistic particles will pass inside the horizon of the black 
hole and be captured; even particles that pass close to the horizon will be captured (Misner, Thorne, 
&: Wheeler 1973; hereafter MTW). As mass accumulates in the black hole, its horizon scale and 
capture radius grow linearly with mass. The pericenter of particles in ballistic orbits, falling from 
our assumed mass distribution, increases as p oc ~ (equation [8] ) . In the earliest stages of 
the collapse, all of the falling material is captured by the black hole. Later, this growth mechanism 
tapers off when the black hole mass reaches a critical point defined by equating the pericenter p 
(for 9 = Tr/2 orbits) to the capture radius of the black hole. In Schwarzschild geometry, particles 
coming inwards from infinity on zero energy orbits are captured by the black hole if p < ARs 
(MTW) , where Rs is the Schwarzschild radius (equation [6] ) . The condition p = ARs defines the 
critical mass scale Mc, 

which represents the mass at which direct accretion onto the black hole becomes compromised. In 
the original version of our model (Paper I), this critical mass scale Mc determined the observed 
black hole mass M^h- Notice that equation [15] displays the correct (observed) scaling with velocity 
dispersion a. 

Even after the centrifugal barrier grows larger than the Schwarzschild radius, however, the 
black hole can gain additional mass from material falling on streamlines that are oriented along the 
rotational poles of the system. After the critical point (described above) is reached, the fraction 
of the infalling material that lands at such small radii is a relatively rapidly decreasing function of 
time. As a result, this effect makes the black hole mass larger by a modest factor J^a, where the 
maximum value !Fa ~ 1-35, as we calculate next: 

The mass infall rate for material falling directly onto the black hole itself is given by 

Mbh= / dn47r{aRsf\vr\p{aRs,n), (16) 
Jo 

where = cos 9. We evaluate the density at the capture surface given by aRs, where a = 8. The 
fact that the capture radius is larger than the Schwarzschild radius is due to the curvature of space 
by the black hole and is a standard relativistic effect (see MTW). Using equation [14] to specify 
the density, we can evaluate the integral to obtain a differential equation for the time evolution of 
the black hole mass, i.e., 

^ = Mbh = M(l - Mc) , (17) 

where fic is the cosine of the angle of the last streamline (measured from the rotational pole) that 
falls directly onto the central black hole. Using the orbit equation [7], we find 

HC = {1- aRs/Rc f^ , (18) 

where we have evaluated /zc at the black hole surface (keep in mind that /xc = for Rc < Rs)- 
Next we define dimensionless variables 



/ = Mbh/Mc and ^ = M/Mc . 



(19) 
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The reduced black hole mass / obeys the ordinary differential equation 

^ = i-{i-f/ef\ (20) 

Equation [20] must be integrated subject to the boundary condition / = 1 at ^ = 1 (which means 
physically that the black hole mass M^h = Mc when the centrifugal barrier Rc first exceeds the 
Schwarzschild radius Rs, which is true by definition). If no additional physical effects prevent 
continued accretion onto the black hole, the enhancement factor is determined by numerically 
integrating equation [20] to find the limiting value as ^ — > oo. In this limit, the black hole mass 
grows by an additional factor Ta as infall continues along the rotational poles of the system, where 

jr4 = lim/(^) 1.3502. (21) 

In other words, the final black hole mass ~ J-aMq, with J-'a ~ 1.35, in the absence of additional 
physical effects. 

However, two additional processes can affect this prediction. First, the infalling material can 
experience shell crossings and the baryonic material can be heated to the system's virial tempera- 
ture. As a result, a hot corona forms around the black hole and the infalling material must cross 
this corona in order to become part of this black hole (anonymous referee; private communication). 
This effect acts to reduce the effective value of J^a and hence equation [21] provides an upper limit. 
In other words, the enhancement factor J^a is confined to the range 1 < J^a < 1-35. Second, 
additional material can be added to the black hole through disk accretion; this process is addressed 
in §3.1. 

We can now evaluate the mass scale for forming black holes using our fiducial value of the 
initial rotation rate (see §2.1), equation [15] to set the critical mass scale, and equation [21] to 
specify the enhancement factor J^a- We thereby find the M^h — cr relation in the form 

Mbh = ^A-^ ~ 1O«M0((7/2OO km s"^)^ , (22) 

where we have written the result in terms of a rather than Oeff . This relation is in reasonably good 
agreement with the observed correlations (see equation [1], Figure 1, and Paper I). 

Scatter in the value of the initial rotation rate $7 will produce corresponding scatter in the 
resulting — cr relation. To obtain an estimate of this effect, suppose that the distribution of CI 
has the same form as the distribution of the spin parameter A for galactic halos, where numerical 
simulations suggest that P(A) oc exp[— (In A/Ao)^/2(7^], with ax ^ 0.5 (Bullock et al. 2001). If the 
initial rotation rate in this model has the same variance, then the ~ c relation will develop 
scatter at the level of (T;^/lnl0 = 0.22 dex. This level of scatter is represented in Figure 1 by the 
bold-faced error bar symbols on the theoretical curve. For completeness, we note that variations 
in the enhancement factor J^a will introduce additional scatter into the Mbh — cr relation. If the 
parameter J^a is uniformly distributed over the range 1 < J^a < 1.35, the resulting scatter will 
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be about 0.04 dex. Since the combined variances of $7 and Ta add in quadrature (e.g., Richtmyer 
1978), the resulting scatter is approximately 0.224 dex (safely smaller than the observed scatter of 
0.30 dex). 

Most of the baryonic material not captured by the black hole during this early collapse phase 
eventually forms stars in the galactic bulge. Dark matter with low angular momentum is captured 
into the black hole along with the baryons; dark matter with high angular momentum {p > 4:Rs) 
passes right through the galactic plane and forms an extended structure. 



2.4. Bulge Mass Scale and Mass Ratios 

This simple dynamical model also predicts a mass scale Mb for the bulge itself: In the ab- 
sence of additional physical processes, the collapse of a structure with initial conditions described 
by equation [2] will produce a "bulge structure" with a well-defined mass scale. Unfortunately, 
bulge formation is complicated by a host of additional processes (see, e.g., Kauffmann, White, 
& Guiderdoni 1993; Cole et al. 1994; Somerville & Primack 1999; and especially Kauffmann &: 
Haehnelt 200). The baryonic gas must cool to form the bulge, and the cooling time can be quite 
long (much longer than the collapse time scales - see the following subsection) . Additional gas can 
be expelled from the bulge through the action of a galactic wind. The dark matter can undergo 
violent relaxation and only some fraction of the dark matter will remain within the bulge itself. 
Finally, a significant fraction of the baryonic material can form the inner portion of the galactic 
disk, rather than remain in the bulge. In spite of these complications, it is interesting to find the 
mass scale Mb defined by the dynamical model alone: 

If the initial protobulge structure is rotating at angular velocity 0, then only material within 
a length scale R = a^s/^ can collapse to form the bulge. Material at larger radii, r > i?, is already 
rotationally supported and will not fall inwards. In the absence of the aforementioned additional 
processes, the length scale R thus defines an effective outer boundary to the collapsing region that 
forms the bulge (see Terebey et al. 1984 for a detailed calculation of how the rotating, collapsing 
inner region can match smoothly onto the static, uncollapsing outer region). The boundary i?, in 
turn, defines a mass scale for the bulge, Mb ^ M{R), i.e., 

Mb = TbJ'dm-^ « 3.3 x lO^°M0^DM(fT/2OO km s-^)^ , (23) 

where the second approximate equality assumes J- dm = 1 and J-b ~ 7r/2. The factor Tb takes 
into account the fact that material along the poles can fall into the bulge even though material in 
the equatorial plane is rotationally supported. The value Tb = 7r/2 is determined by integrating 
the initial density distribution (equation [2]) over the entire cylinder defined hy w < R = a^s/^ 
(where w is the cylindrical radius). The fraction Tdm < 1 is the fraction of the initial mass that 
is retained within the bulge structure; for example, not all of the dark matter necessarily remains 
in the bulge. 
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The resulting expression (equation [23]) shows that the bulge mass scale exhibits power-law 
behavior with Mb ~ ct^. It is interesting to compare this result to the observed bulge masses in 
the sample of galaxies with central black holes. Figure 2 shows the mass scale of equation [23] 
along with the observed data; also shown is the best fit power-law, which has slope ~ 3.3 ±0.1. For 
comparison, the traditional Faber-Jackson relation for elliptical galaxies and bulges has the form 
L oa (Binney &; Merrifield 1998; Faber & Jackson 1976). Consistency of equation [23] with the 
Faber-Jackson law would require that L oc M^^^. For a constant stellar IMF, this relation, in turn, 

1/3 

would imply a star formation efficiency e of the form e ~ . In other words, in order for the 
mass scale of equation [23] to describe the masses of observed bulges, the more massive systems 
would have to be more effective at forming stars. 

We are also interested in the ratio of black hole mass to bulge mass. In our simple dynamical 
model, the bulge mass scale Mb and the black hole mass Mbh ~ Mc have the same functional 
dependence on the rotation rate The resulting ratio Mbh/M^ = Ms is independent of and is 
given by 

,j^B = ^ = ^-^ 0.0024 (^7/200 km s"^) . (24) 
Mb j^b c 

This mass fraction ^b is roughly comparable to the observed ratio of black hole masses to bulge 
masses in host galaxies. The first estimates suggested that this mass ratio is nearly constant 
(e.g., Richstone et al. 1998; Magorrian et al. 1998), although the data show appreciable scatter. 
Later papers found values of /is = 0.0015 - 0.0020 (e.g.. Ho 1999; Kormendy 2000), in reasonable 
agreement with the typical value suggested by equation [24]. More recent work (Laor 2001) finds 
that the mass ratio is an increasing function of bulge mass, jiB oc M^'^^^*^'^'^, with jib ~ 0.0005 for 
the smallest bulges with observed black holes and hb ^ 0.005 for bright ellipticals (cf. McLure & 
Dunlop 2002); this latter result is somewhat steeper than the prediction of this model, which implies 
//B oc M^"^ . The law [24] is shown in Figure 3 along with the observational data. An unweighted 
fit to the data implies a slope of ~ 0.9 (close to the model prediction of 1.0), but the error bars and 
scatter in the data are too large to make a definitive claim. Notice also that the observed black 
hole scaling law, Mbh ~ cr^, and the observed scaling law for bulge masses, Mb ^ <7^'^, imply that 

7 

fiB ^ a"-. 



2.5. Time Scales 

With an initial density profile of the form p r^^, a detailed collapse calculation (Shu 1977) 
indicates that the flow exhibits a well defined mass infall rate M = moa^g/G, where mo « 0.975 
[notice that this starting state corresponds to an unstable hydrostatic equilibrium] . The infall rate 
is constant in time and we can measure the time elapsed since the collapse began by the total mass 
M that has fallen to the galactic center. At early times, all of the mass falling to the center is 
incorporated into the central black hole. At later times, the mass supply is abruptly shut off by 
conservation of angular momentum. In this setting, the mass infall rate is quite large, M 650 
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Mq yr~^ (for a = 200 km s~-'^and = Sa^g). The time scale Tbh to form a typical supermassive 
black hole (with mass ~ 10^ ^'^o) is about Tbh ~ 10^ yr, comparable to the time scale r* 

for individual stars to form (e.g., Adams & Fatuzzo 1996; Myers & Fuller 1993). In the absence 
of any competing physical processes, the dynamical time scale to form the entire bulge is much 
longer, about Tbig ~ 25 — 50 Myr, comparable to the crossing time tcross = R/c^eS of the initial 
structure. One should keep in mind that bulge formation also requires the gas to cool, however, 
and the cooling time scale can be longer than this dynamical time scale. 

3. GENERALIZATIONS OF THE MODEL 

In this section wc present further generalizations of this collapse model for black hole forma- 
tion. The previous section shows how the black hole gains mass through infall from the collapse 
flow. However, additional mass can be added to the black hole through disk accretion (§3.1). Fur- 
thermore, black holes formed through this collapse picture are born rapidly rotating (§3.2). The 
scaling relation for the black hole mass as a function of a depends on the initial angular 
momentum profile of the pre-collapse material; in fact, the angular momentum distribution is the 
most important determining factor in specifying the black hole masses (see §3.3, 3.4). This model 
can also be generalized include the effects of mergers (§3.5) and non-spherical initial conditions (in 
particular, quadrapole moments; see §3.6). 

3.1. Disk Accretion 

Material that falls to the niidplanc of the system in gaseous form can collect into a disk 
structure that surrounds the nascent black hole. The presence of the disk is consistent with the 
current theoretical ideas about active galactic nuclei and the jets they produce. In order to retain 
the desired scaling law Mbh ~ cr^, however, the total amount of mass added to the black hole 
through disk accretion should be less than (or at most comparable to) the original mass scale Mc- 

The energy density of the universe from quasar activity places a limit on the amount of 

mass that can be accreted by black holes. This energy density Ut has been estimated to be 
Ut ~ 2.5 X 10^^^ erg/cm^ (Elvis, Risaliti, & Zamorani 2002). If this energy arises from mass 
accretion onto black holes with energy conversion efficiency e, then the energy density Ut implies 
a corresponding minimum mass density in black holes at the present cosmological epoch. Two 
recent papers (Elvis et al. 2002; Yu &; Tremaine 2002) have estimated the amount of mass accreted 
through quasar activity over the reshift range < z < 5. The mass density in black holes due to 
accretion activity takes the form p^cc ~ 2.1 x 10^ [0-1(1 — Mq Mpc~^. For comparison, the 
observed Mbh — f relation implies a present day mass density in black holes pbh ~ 2.5 x IO^Mq 
Mpc~^ (Yu & Tremaine 2002; Aller & Richstone 2002). The present ratio of the accreted mass to 
the observed mass is thus TZ = 0.083 (1 — e)/e. In order for the infall collapse model of this paper 
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to explain the observed black hole masses, the ratio TZ < 1 and hence the conversion efficiency 
must be at least e ~ 0.08. If the efficiency of energy conversion in quasars is too low, then black 
holes would gain more mass from disk accretion than from infall. If the energy conversion efficiency 
is high, however, quasar light could be explained by a modest addition of mass, and some other 
physical process (e.g., infall) would be required to explain the observed black hole masses. 

Because of the present observational uncertainties, the critical value of ec (the value required 
for disk accretion to explain the observed mass density in black holes) is not completely specified. 
Using only the quasar constraint, Yu & Tremaine (2002) find ec 0.1, whereas Elvis et al. (2002) 
find ec ~ 0.15. For Schwarzschild black holes, the energy conversion efficiency is expected to be 
about e ~ 0.1, but higher efficiency, with e ~ 0.2, is possible for thin-disk accretion onto a Kerr 
black hole (see Yu & Treaminc 2002). A plausible upper limit for accretion processes is e = 0.31 
(Thorne 1974). Both the observations of background radiation and the theoretical expectations for 
energy conversion efficiency e should be specified further to resolve this issue. 

For comparison, we can derive another constraint on the amount of accreted mass. This 
constraint is more general than in the discussion above, but is also weaker. We consider the limiting 
case in which disk accretion is maximally effective. Disk accretion generally cannot operate faster 
than the orbit time at the outer disk edge (Shu 1992). In this context, the orbit time r is given by 

^ -~GM~- ■ ^^^^ 

When the disk accretion time becomes longer than the time required for the disk to condense into 
stars, i.e., when r > r*, disk accretion is no longer effective and the mass flow onto the central black 
hole must come to an end. This condition implies a maximum mass scale for accreting black holes. 
The mass appearing in equation [25] above represents the total mass that has fallen to the center 
by a given time. Only a fraction of this mass is available to join the disk because only a fraction 
fs is baryonic and a fraction fc is in gas (rather than in stars). Including these two factors, the 
maximum mass that can be added to the black hole via disk accretion becomes 

M^ax = (27r)-VVB/G(c7VG)(r,/J^3)i/4 . (26) 

For typical values of the input parameters o and Jl, and for = 2/15, fo = 1/2, this maximum 
mass scale is about 5 times larger than Mc- In the limit of maximally efficient disk accretion, 
the mass scale Mq defined by the centrifugal radius can thus be compromised (see also Burkert &; 
Silk 2001; Silk h Rees 1998). Unfortunately, the disk accretion rates are not known in these early 
stages of galaxy formation (however, see Kumar 1999). In most known astrophysical disks, the disk 
accretion rates are much smaller than their maximum values, typically by factors of 10^ — 10^ (e.g., 
Shu 1992), which would imply that most of the black hole mass does not come from accretion. 
On the other hand, as discussed above, the observed X-ray background implies that the central 
black holes that power quasars must accrete a substantial mass, roughly comparable to the masses 
obtained via infall. This argument assumes that the mass in the disk is large enough to affect the 
black hole mass; however, the infall model used here naturally places most of the incoming mass 
at large radii appropriate for the disk (e.g., Cassen h Moosman 1981). 
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3.2. Black Hole Angular Momentum 



This theoretical model for supermassive black hole formation predicts that the resulting black 
holes should be rotating rather rapidly. If we approximate the orbital solutions for incoming 
material using the classical treatment of §2.2 and the effective capture radius of a Schwarzschild 
black hole, then the angular momentum Jbh of the resulting black hole takes the form 



In the absence of additional mass sources, we can write this angular momentum in terms of the 
black hole mass Mbh = Mc and the Schwarzschild radius Rs, 



The maximum allowed value for the numerical coefficient in equation [28] is 0.5 (e.g., Thorne, Price, 
& MacDonaold 1986; Blandford 1990), so these black holes are rotating close to their maximum 
rates. Notice that relativists usually write the Kerr metric in terms of the parameter a = J/M. 
In units where G = 1, this parameter has a maximum value of a = M. This theory predicts 
the formation of black holes with the ratio a/M 0.89 - close to its maximum value of unity. 
Continued infall will add both mass and angular momentum to the black hole, and will change this 
prediction somewhat. Notice also that disk accretion models predict that supermassive black holes 
should be rotating rapidly (e.g., Elvis et al. 2002). 

Keep in mind that this result uses mixed approximations. We have derived the maximum 
mass scale using the capture radius appropriate for Schwarzschild geometry and purely classical 
orbits solutions. The resulting black holes are rapidly rotating, however, and curve space-time as 
described by the Kerr metric. A fully self-consistent calculation should find orbit solutions and 
the capture radius for Kerr black holes, with the angular momentum parameter aj = Jbh/^^bh 
determined at each evolutionary stage by conservation of angular momentum. This calculation is 
beyond the scope of this present paper, but could change our results at the level of 50 percent. 

The last stable orbit for the Kerr metric can be worked out in terms of simple functions (e.g., 
Riffert 2000) ; the last stable orbit for Kerr geometry is smaller than that of Schwarzschild geometry, 
and the last captured streamline should behave similarly. In other words, a Kerr black hole has a 
smaller surface area than a Schwarzschild black hole of the same mass (e.g., Rees 1984). As the 
black hole grows, it gains both mass - which makes its capture cross section larger - and angular 
momentum - which makes its capture cross section smaller. Near the end of the infall phase, the 
added angular momentum to the black hole decreases its cross section, and the infalling material 
starts to have too much angular momentum to fall within the older, larger cross section. This effect 
thus acts to make the transition more abrupt. 




(27) 



Jbh = ^cMi^hRs « O.UcMbhRs ■ 



(28) 
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3.3. General Initial Conditions 



The success of this simple model for black hole formation rests on the initial angular momentum 
profile, which must have a particular form. In order to determine how sensitive the results are to 
the assumed initial conditions, we consider generalized initial density distributions of the form 



p(r) = Ar 



(29) 



where ^ is a constant and F is an arbitrary index. As before, we assume that the initial structure 
is rotating at a uniform rate J7. The corresponding mass distribution is given by 



M(r) 



4ttA 
3^' 



Ar 



3-r 



(30) 



where we have defined a reduced constant A = An A/ {3 — F). To find the black hole mass scale 
resulting from the collapse of this initial state, we use the same criterion as before, which takes the 
form joo = 4GM/c, where joo = ^ ^lo- To specify the starting radius roo, we must invert the mass 
distribution to obtain 

r Ml l/f3-r) 

(31) 



A 



Solving for the black hole mass scale, we obtain 



M, 



bh 



4G 



(3-r)/(r-i) 



^2/(r-i) ^ 



(32) 



where we must restrict the analysis to F > 1. For less steep initial density distributions, the 
centrifugal barrier of the collapse flow increases more slowly with mass than does the capture 
radius of the black hole (which has a linear dependence) . 

We can also solve for the other parameters of the forming bulge system. The radius R that 
marks the outer boundary is given by 



R 



AG 



i/r 



while the bulge mass scale takes the form 

Mb -- 



G 



Ln2 



3/r-i 



(33) 



(34) 



With the radius R and bulge mass scale defined, we can solve for the expected velocity dispersion 
of the final system using the relation 2GMb / R, where the factor of 2 arises from the collapse 
itself. The resulting expression for the velocity dispersion is 



a = V20i-2/r 



AG 



i/r 



(35) 
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With this expression for a in hand, we can write the expected bulge mass Mb, bulge size scale R, 
and black hole mass M^h in terms of a rather than the initial variable A (or A) appearing in the 
density distribution, i.e.. 



.3 ^4 



(4-2r)/(r-i) 



Mb = , R = a/{V2n) , and M^h = ^ \V8- " " " " . (36) 



And finally, we obtain the corresponding expression for the mass ratio: 



V8- 



(3-r)/(r-i) 



(37) 



This result is very sensitive to the starting slope F of the density profile (for an assumed constant 
rotation rate Q). If the index V is much smaller (larger) than the preferred value T = 2, then the 
final black hole masses are much smaller (larger) than those observed. Variations in the slope F 
will produce corresponding scatter in the observed — a relation; as a benchmark, variations at 
the level of F = 2 lb 0.1, will produce scatter of about 0.5 dex. The value of the index F affects this 
model for mass scales Ms in the range O.lMbh < Ms < 0.5Mb. 

In summary, this model is sensitive to the form of the initial conditions. Moderate departures 
from our assumed starting condition p ~ to large variations in the predicted black hole masses 
-Mbh and mass ratios iib- If this model is correct, then the relevant initial conditions for the collapse 
of galactic bulges must have angular momentum distributions close to those assumed here (specified 
by equation [2]). 



3.4. Dimensional Analysis 

Given all of the generalized relations discussed in the previous section, how can we make sense 
of all the possibilities? For this theoretical model, the story is relatively simple: For the formation 

of the bulge itself, the black hole forming at the galactic center has essentially no effect. Because 
no relativistic effects come into play, the speed of light c does not enter into the formulae describing 
the bulge properties u, Mb, and R. The only variables that can determine these quantities are the 
rotation rate Vl, the parameters of the initial density distribution [A and F), and the gravitational 
constant G. Furthermore, the index F is dimensionless, so the only quantities that carry dimensions 
arc A, n, and G. A simple virial argument lets us replace the density coefficient A with the velocity 
dispersion fj, which is a much more familiar quantity. For a given velocity scale cj, the quantity 
a'^ /G defines a mass infall rate; when combined with a time scale (which is the only time scale 
present in this simple treatment) we thus obtain the mass scale M ~ a^/GQ, which we identify 
as the bulge mass scale Mb- To summarize, the three quantities with dimensions {A, G, $1) can 
only define the three bulge quantities (cr, Mb, R) in one way (up to dimensionless factors of order 
unity) . 

Next, however, we must put the black hole mass into the problem. The black hole intro- 
duces relativistic effects and, in particular, introduces the speed of light c as another fundamental 
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constant. The ratio (3 = a/c is another dimensionless parameter in the problem. In terms of di- 
mensional analysis, we can now define an infinite number of new scales. For example, in addition to 
the bulge mass scale Mb, we now have the family of mass scales M = F{P)Mb, where F{P) is an 
unspecified function of the second dimensionless field (3. In our treatment, the power-law forms for 
the initial conditions lead to new mass scales given by the power-laws = Mb P^, where rj can 
be any real number. The physical law of angular momentum conservation, in conjunction with the 
initial profile, chooses the value of the exponent f] for a given model. In the simplest case, that of 
isothermal p r^^ initial conditions, we obtain rj = 1 and hence Mbh ~ Mb{(t/c), which happens 
to be the observed scaling law (Mbh ~ f'"^)- Notice also that the ratio Afbh/M^ ~ tj/c ~ 10~^ (the 
observed order of magnitude). The other (more general) choices of initial conditions, p ~ r^^, thus 
correspond to other possible values of the exponent rj, i.e., = (3 — r)/(r — 1). 



Many galaxies are expected to experience merger events during their formative stages of evo- 
lution (e.g.. White & Rees 1978; White 1979). If the initial collapse of protobulges proceeds as 
envisioned in the simple theoretical model developed in this paper, then the fundamental building 
blocks of galaxies will have bulge mass Mb ~ o"'^ and black hole masses Mbh o"*. However, 
observations indicate that the final merger products obey these scaling relations. As a result, we 
must now consider what happens to these scaling laws if the fundamental building blocks undergo 
multiple merger events (see also Hughes & Blandford 2002; Menou, Haiman, &; Narayanan 2001; 
Ebisuzaki et al. 2001; Volonteri, Haardt, & Madau 2003). 

For the sake of this discussion, we consider the rotation rate to be constant so that the only rel- 
evant variable is the velocity dispersion a. We define a dimensionless variable s = cr/(200 km s~^). 
If protobulges merge many times and if they all contain central black holes that merge, then the 
final black hole mass can be written in the form 



where we assume in the second equality that the starting units, labeled by the index j, obey the 
scaling law derived in §2 (and Paper I). All sums are taken up to N, which specifies the number 
of protobulges (initial units) that merge to form the final stellar system. Notice that this equation 
assumes minimal energy losses from gravitational radiation during the collisions. 

We must relate the velocity dispersion s/ of the final bulge system to the velocity dispersions 
Sj of the individual units. We first consider one idealized limiting case: If the mergers take place 
with zero orbital energy and no energy losses occur during the collision, then the final energy of 
a merged system must equal the internal energies of the initial (pre-merger) pairs. If we consider 
the systems to be in virial equilibrium and to be homologous (a gross approximation, but a good 
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(38) 



J 



J 



-19- 



place to start for conceptual purposes), we obtain the relation 



J2pjs], (39) 



where we have defined pj = RBj/Rsf- 



Now let us define qj = s| for all protobulges labeled by the index j. We can think of the qj and 
the Pj as components of A'^-dimensional vectors (which live in the space of protobulge parameters). 
With this ansatz, we can write the final black hole mass in the form 

Mbh f = Mo^q] = Mo\qj\\ (40) 
j 

where the notation \qj\'^ denotes the vector magnitude squared; the final velocity dispersion Sf takes 
the form 

s} = p q. sj = \pj\'^\qj\'^ cos'^ 01,1, , (41) 

where ^?bh is the angle between the two vectors in bulge parameter space. Combining the above 
equations, the final black hole mass takes the suggestive form 

The final mass of the black hole thus displays a quartic dependence on the final velocity dispersion 
(to leading order). Some intrinsic scatter will be introduced through varying angles ^^h ^-nd mass 
weights Pj. We can consider the limiting case in which all of the starting protobulge units are 
identical, pj = 1/N for all units j, and the two vectors pj and qj are parallel so that cos^bh = 1- 
In this limit, the final black hole mass takes the form Mbh = NMqs^, which is consistent with the 
observed scaling relation (equation [1]). In other words, the scaling law can be preserved under the 
action of mergers in this idealized limit. 

We can also find the behavior of the bulge mass under repeated merger events. The final bulge 
mass Msf after N mergers is given by 

Mb f = MBoY^ 4 = ^jlj = ^Bo s • q . (43) 

3 3 

We can write the dot product in terms of another angle Ob, 

s • q = \qj\ cos 9b , (44) 
which allows us to write the final bulge mass in the form 
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To explicitly illustrate the possibility of preserving the M^h — cr relation under the action of 
mergers, we have performed a simple set of Monte Carlo simulations, as depicted in Figure 4. In 
these numerical experiments, the number of interacting units is randomly chosen within the range 
2 < iV < 6 (see, e.g., Wechsler et al. 2002). The velocity dispersions of the interacting units are 
chosen to be random, but logarithmically spaced in order to evenly sample the observed range of 
(final) a. The radial sizes of the interacting units = RjB/RfB are chosen to be equal so that 
Pj = 1/N. In order to get the normalization correct, we let the mass scale Mq of the interacting 
units (equation [38]) be a factor of three smaller than that of the observationally determined 
distribution (equation [1]). Notice that this approximation is equivalent to allowing serious energy 
losses during the merger events. The final values of the black hole mass and velocity dispersion are 
then calculated according to equations [38] and [39]. With these idealizations, the final distribution 
retains its power-law form and appears to be in good agreement with the observed correlation (see 
Figure 4). We stress, however, that the range of parameter space available for merging protobulges 
is enormous - not every merger scenario will produce such a clean correlation. We leave a more 
detailed exploration of parameter space for future work. 



3.6. The Effects of Initial Quadrapole Moments 

In this subsection, wc consider the effects of a quadrapole moment on this collapse picture of 
black hole formation. Since the protobulge structures can obtain their initial rotation rates through 
tidal torques acting on nonzero quadrapole moments of the mass distribution, it is important to 
check whether such quadrapole moments affect the collapse flow. We can consider two limiting 
cases: A substantial quadrapole moment in the inner regime and a substantial quadrapole on the 
large size scale of the protobulge itself. The case of the inner regime is very much like the binary 
star potential in the star formation problem, and this situation has already been shown to have 
little effect (Allen 1999, Jijina 1999). We thus consider the case of an outside quadrapole moment. 

As a starting point, wc consider the outer quadrapole to be akin to two point masses M with 
separation R. In order for the outer quadrapole to exert a torque on the inner region (the portion 
of the pre-coUapse structure that will eventually comprise the black hole), the inner region must 
also have a quadrapole moment. We consider the inner portion to have a dumbbell shape with 
mass scale m and size £. The force exerted on the inner region by the outer quadrapole is given by 
the tidal force law 

^ GMm e 

The torque r exerted on the inner region is given by this force acting through a lever arm of size 
£. The torque is given by 

r-^^. (47) 

Now we need to make the connection between these quantities and the scales in the black hole 
formation problem. The size scale R is the size of the protobulge, i.e., R Qes/^- Because the 
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mass distribution has the profile of an isothermal sphere {p ~ and M(r) ~ r), the size i is 
smaller than R by the ratio of the black hole mass to the bulge mass, i.e., d/R = fis- The mass M 
of the outer quadrapole is given by M = XiMb, where the dimcnsionless parameter Ai determines 
the asphericity of the configuration. If the distribution is perfectly spherical, Ai 0. In the limit 
that the protobulge has a dumbbell shape, Ai —>■ 1. Similarly, we let m = A2Afbh- With these 
definitions, the torque can be written in terms of the physical quantities in our paper as 

T = \1\2HB ^ • (48) 

One way to quantify the effect of this torque on the collapse flow is to find the total change 
in angular momentum produced by the torque, which acts over the collapse time At of the inner 
region only. (After the black hole has formed, the still-existing outer quadrapole will not effect the 
black hole mass). The time for the inner region to collapse is given by 

Combining the above expressions for the torque and the time interval over which it acts, we 
find the change in angular momentum: 

AJ = 2AiA2/.|^^ . (50) 

For later convenience, we note that our model implies GMb/R^ = 2CP, so we rewrite the above 
expression to obtain the form 

AJ = 4AiA2AiBMbhJ^^^ . (51) 

The starting angular momentum Jq of the inner region can be obtained by integrating over 
the initial density distribution. The result is 

Jo = ^Mki^nf . (52) 

The resulting fractional change in the angular momentum of the inner result is immediately found 
to be 

AJ 

— = 18AiA2/xs. (53) 

The mass fraction ijlb ^ 0.0024 in our model, so this fraction change becomes 

AJ 

— «i O.O43A1A2 . (54) 

Even in the extreme limit of highly developed quadrapole moments on both the inside and the 
outside, the fractional change in angular momentum is only about 4 percent. In a more realistic 
case, the outer quadrapole should be smaller than unity Ai < 1, but still large enough to give the 
bulge an elliptical appearance; the inner region can be smoothed out (e.g., see Peebles 1993 for 
streaming arguments) and the inner quadrapole should also have A2 < 1. In any case, the coupling 
between the outer quadrapole moment and the inner collapse region can be considered weak in the 
context of our orbit solutions. 
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4. CONCLUSIONS 

In this contribution, we have presented further development of the simple theoretical model 
for supermassive black hole formation that was put forth in Paper I. This model begins with an 
initial state specified by a density distribution of the form p = a^^/2iTGr'^ and a uniform rotation 
rate The parameters (aefr,ri) represent the specification of the initial conditions. As the initial 
state collapses to form a galactic bulge, the collapse flow produces a black hole in the center. The 
velocity dispersion of the final bulge system is comparable to the effective transport speed and we 
make the identification a ~ \/2fflcfr- In developing this basic picture, we find the following results: 

[1] The black hole mass Mbh is determined by the condition that the centrifugal radius ex- 
ceeds the capture radius of the central black hole. This requirement leads to the scaling law 
Mbh = Mo(c7/200 km s""^)^, which is consistent with observations both in its dependence on veloc- 
ity dispersion a and for the mass scale (Mq ~ IO^Mq) of the leading coefficient (see equation [22] 
and Figure 1). The mass scale Mq is given by Mq = 4.Fyi(200km/s)^/cGr2, so that variations in 
the rotation rate and the amount of continued infall J- a lead to dispersion about the observed 
power-law correlation. If the initial rotation rate follows the same distribution as that calculated 
for the spin parameter A of dark matter halos, then variations in Q. would produce a scatter of 
~ 0.22 dex (a factor of ~ 1.7) in the predicted Mbh — cr scaling law. The observed relation has a 
factor of 2 dispersion. 

[2] A bulge mass scale is defined in this model by the outer boundary of the collapsing region. 
Material at initial radii r > R is rotationally supported and cannot collapse, where R = a^s/^- 
This condition implies a scaling law for the bulge mass scale, Mb = 2!FBal{[/GCl oc (see equation 
[23] and Figure 2). Although bulge formation must involve physical processes that are not included 
in this dynamical mdoel (e.g., gas cooling, feedback from galactic winds, disk formation), this 
scaling law for Mb is nonetheless in good agreement with the observed relation for host galaxies 
that contain supermassive black holes, as shown in Figure 2. 

[3] If we interpret the bulge mass scale Mb (see [2] above) as the bulge mass itself, then this 
model predicts the ratio p,B oi black hole mass to bulge mass (equation [24] and Figure 3). The 
theoretical mass ratio is proportional to the velocity dispersion and has the form hb 0.0024 
(cr/200 km s~^), roughly comparable to observed mass ratios. A mass ratio }ib that increases with 
velocity dispersion a is consistent with (and even indicated by) a recent analysis of the observational 
data (Laor 2001); an unweighted least squares fit to the data shown in Figure 3 implies a slope of 
~ 0.9. We also note that a constant mass ratio jiB may be inconsistent with the data: Since the 
black hole mass scales as M^h ~ cr^ (Tremaine et al. 2002), a constant mass ratio iib would require 
the bulge mass to also scale as Mb ~ cr^. However, the observed bulge masses (for systems with 
detected black holes) do not indicate such a steep dependence on cr; the data suggest an index of 
approximately 3.3 it 0.1 (see Figure 2). 

[4] In this scenario, the black hole forms quickly, with a typical formation time of ~ 10^ yr. 



-23- 



[5] The black holes formed through this mechanism are born with rapid rotation rates. As 

a result, the geometry in the central regions is best described by the Kerr metric. Specifically, 
this model predicts that supermassive black holes arc formed with an initial rotation parameter 
a/M K, 0.9, relatively close to the maximum value of a/M = 1. Such high black hole rotation rates 
may be detectable by LISA (the Laser Interferometer Space Antenna). 

[6] Although the supermassive black holes produced by this process are intrinsically relativistic 
objects, relativistic effects play only a modest role in the collapse flows that produce them. The 
most important effect is that the capture radius of a black hole is larger than the Schwarzschild 
radius by a factor of 4 and this effect leads to black hole masses that are larger by this same factor. 

[7] The predicted black hole masses are very sensitive to the initial conditions. For an initial 
density distribution of the form p ~ (with constant angular velocity O), the subsequent collapse 
produces black holes and galactic bulges with the correct masses and the correct dependence on 
velocity dispersion a. The correct mass normalization depends on the choice of rotation rate 
f2 ~ 2 X 10~^^ rad/s. Initial density profiles with shallower slopes, p ~ with F < 2, produce 
smaller black holes with a steeper slope in the log Mbh — log a plane. In general, the initial profile of 
specific angular momentum - given by the combination of p{r) and ^{r) - determines the final mass 
scales. This sensitivity on initial conditions is both the strongest and weakest aspect of the model: 
If we can unambiguously determine the angular momentum distribution of the initial states, we 
can directly verify or falsify this theoretical scenario for black hole and bulge formation. We also 
note that our initial conditions apply on (initial) radial scales of several kpc. The manner in which 
these initial conditions match onto the Hubble flow and the larger scale structure of the galactic 
halo (at radial scales of many hundred kpc) remain to be determined. 

[8] If galactic bulges and larger galactic structures are formed through the mergers of smaller 
constituent pieces, this scenario for black hole formation can still play a role: In this case, a num- 
ber of the constituent pieces would form black holes in their centers through this mechanism. The 
resulting scaling laws (e.g., Mbh ~ cr*^) can be preserved, or nearly preserved, under mergers for 
idealized circumstances (§3.5). On the other hand, mergers tend to reduce the angular momen- 
tum per unit mass, so that merger scenarios predict a lower angular momentum for the resulting 
black holes (e.g., Hughes & Blandford 2002). Future measurements of the angular momentum of 
supermassive black holes are thus crucial for discriminating between various formation scenarios. 

We would like to thank Gus Evrard, Karl Gebhardt, and Risa Wechsler for useful discussions; 
we especially thank Rosie Wyse for enlightening discussions regarding bulge rotation rates. Finally, 
we thank an anonymous referee for many suggestions that improved the paper. This work was 
supported by NASA through the Long Term Space Astrophysics program and the Space Telescope 
Science Institute, and by the University of Michigan through the Michigan Center for Theoretical 
Physics. 
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Fig. 1. — The correlation between blaek hole mass M^h and velocity dispersion a of the host 
galaxy. The data points (adapted from Gebhardt et al. 2003) represent the observed correlation 
for ellipticals (circles), SO galaxies (squares), and spirals (triangles). The solid curve shows the 
theoretical result of this paper (using equation [22] with Ta = 1-35). The dashed curves show the 
observational fit advocated by Tremaine et al. (2002) ; curves are shown for the best estimate of the 
index 7 = 4.02 and for the maximum/minmum values 7 = 4.02 it 0.32. The bold-faced error bar 
symbols show the level of scatter that would result if the initial rotation rate Q follows the same 
distribution as that of the spin parameter A for galactic halos. 
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Fig. 2. — The correlation between bulge mass scale Mb and velocity dispersion a of the host 
galaxy. The data points (adapted from Gebhardt et al. 2003) represent the observed correlation 
for ellipticals (circles), SO galaxies (squares), and spirals (triangles). The error bars correspond to 
20 percent uncertainties in the bulge mass estimates. The solid curve shows the theoretical mass 
scale predicted by the infall-coUapse model of this paper (using equation [23] with JTg = 7r/2 and 
^DM = !)• The dashed curve shows an unweighted least squares fit to the data. 
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Fig. 3. — The ratio /x_b of black hole mass Mbh to bulge mass scale Mg plotted as a function 
of the velocity dispersion a of the host galaxy. The solid curve shows the prediction of equation 
[24], where we assume that the simple dynamical mass scale Mb from the collapse model can be 
identified with the bulge mass. The data points (adapted from Gebhardt et al. 2003) exhibit 
considerable scatter, but the least squares fit (shown as by the dashed curve with slope 0.86) is 
in reasonable agreement with theoretical expectations; the various symbols represent ellipticals 
(circles), SO galaxies (squares), and spirals (triangles). The error bars are determined by the 
quadrature sum of the error bars shown in Figures 1 and 2. 
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Fig. 4. — An illustration of the correlation between black hole mass Mbh and velocity dispersion 
G being preserved under the action of mergers. The result of each numerical experiment is shown 
as an open square. The lines (for reference) are the same as those in Figure 1. The numerical 
experiments begin with N = 2 — ^ smaller units, which are merged according to the rules of §3.5. 



